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Abstract. In this paper we study the integral of the supremum pro- 
cess of standard Brownian motion. We present an explicit formula for 
the moments of the integral (or area) A{T), covered by the process in 
the time interval [0, T]. The Laplace transform of A{T) follows as a con- 
sequence. The main proof involves a double Laplace transform of A{T) 
and is based on excursion theory and local time for Brownian motion. 



1. Introduction 

Let B(t), t > 0, be a standard Brownian motion. Consider the following 
associated processes: the supremum process S{t) := maxo<s<t B{t), and the 
local time L{t), which can be regarded as a measure of the time B(t) spends 
at in the interval [0,t], see Revuz and Yor fid, Chapter VI] for details. It 
is well-known that these two processes, although pathwise quite different, 
have the same distribution [10|, Chapter VI. 2], 

m}t>o={mt>o- 

The purpose of this paper is to study the distribution of the area under 
S{t) or, equivalently, L{t) over a given time interval [0,T]. That is, the 
integral 

A{T):= r S{t)dt= r L{t)dt. (1.1) 
JO Jo 

For ease of notation, let ^ := ^(1). 

The area (ll.ip appeared as a random parameter when analysing displace- 
ments for linear probing hashing. The Laplace transform of A, which is 
presented in Corollary 12.41 provided the means to prove one of the main 
theorems in Petersson [9]. 

Note that the usual Brownian scaling 



{BiTt)}^^^^{T'/'B{t)} 



t>0' 
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for any T > 0, implies the corresponding scaling for the supremum process, 
Thus, for T > 0, 

A{T) = T [ S{Tt) dt = T^/^A, (1.2) 
Jo 



and it is enough to study A. 



2. Results 



Let ip{s) := Ee~*"^ denote the Laplace transform of A. An essential part 
of this paper is devoted to proving the following formula for the Laplace 
transform of a variation of tp, or in other words, a double Laplace transform 
of A. Such formulas have already been derived for the integral of |-B(i)| and 
other similar integrals of processes related to Brownian motion, see Perman 
and Wellner and the survey by Janson [3']. 

Theorem 2.1. Let tp be the Laplace transform of A. For all a, A > 0, 

Remark 2.2. One of the parameters a and A in Theorem 12.11 can be elim- 
inated (by setting it equal to 1, for instance) without loss of generality. In 
fact, for any /? > 0, the formula is preserved by the substitutions A ^ j3\, 
a I— > 0^/'^a and s i— > [3~^s. 

The proof is given in Section [3 It is based on excursion theory for Brow- 
nian motion and is inspired by similar arguments for other Brownian areas, 
see Perman and Wellner [§]. 

Theorem 2.3. The n:th moment of A is 

EA-= "'r(n + 2/3) (3V2y 

r(2/3)r(3n/2 + l)l 4 ) ' ^ 

Proof. Set A = 1 in Theorem 12.11 and denote the left and right hand side by 

/•oo 

1(a) := / ip{as^/^)e-' ds 
Jo 

and 



T/ N ; / -, 3q;s \ -2/3 







The integrand of I{a) and all its derivatives with respect to a are dom- 
inated by functions of the form s^e~*, uniformly in a > 0. Differentiation 
of /(a) is therefore allowed indefinitely due to dominated convergence. The 
same argument applies to J(a). 
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Also, the dominated convergence theorem shows that integration (with 
respect to s) can be interchanged with taking the hmit a 0+. Thus 

d"/(a) ^. p d" 3/2, 
hm — --^ = hm / - — ■ip(as'^'^)e Ms 
a^o+ do" a-*0+ Jq da"^ ^ 

= hrn^(-s3/2)«E(^"exp{-as3/2^})e-Ms 

/•oo 

= (-l)"E(^") / s3n/2g-.^^ 

Jo 



(-l)"r(3n/2 + 1)E^" 



and 



hm — = hm / 1 H e'" ds 

a^o+ da"' a->o+ da" V 2^2^ 



7o r(2/3) V2^/ V 2^2^ 

_ r(n + 2/3) / -3 \" /■ 
r(2/3) 12^2^ Jo 
^ r(n + 2/3) / -SVl y , 

r(2/3) V 4 ; 

The fact that /(a) = J(a) completes the proof. □ 

The first four moments of A are listed in Table [TJ Further, Stirling's 
formula provides the asymptotic relation 

Corollary 2.4. The Laplace transform of A is 

1 r(n + 2/3) f -3V2s Y 

= fm r(3n/2 + 1) l^^J • ^'-'^ 

Proof. The corollary follows from the identity 



n! 

n=0 



Note that the sum converges absolutely for every complex s. □ 

The graph of ip{s) is shown in Figure [TJ 

Remark 2.5. The Laplace transform of A can also be expressed in terms 
of generalized hyper geometric functions, 

/5 4 s2\ 4s /6 8 7 9 s^ 
V'(s) - 1^1 6= "6 J " V2^ U' 6' 6' 6' ¥ 
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3V27r 



=12 



64 



63V27r 



E^^ = — 
24 



Table 1. The first four moments of A. 



1.0 -, 



0.8 - 




Figure 1 . The Laplace transform of A. 



3. Tail asymptotics 

Tauberian theorems by Davies [1] and Kasahara [2] (see Janson [4, Theo- 
rem 4.5] for a convenient version) show that the moment asymptotics (|2.ip 
implies the estimate lnP(^ > x) ~ —3x^/2 for the tail of the distribution 
function. Thus, the following corollary is obtained. 

Corollary 3.1. A has the tail estimate 

> x) =exp{-3xV2 + o(x^)}, x ^ oo. 

(This result can also be proved by large deviation theory; cf. similar results 
in Fill and Janson 

It seems difficult to obtain more precise tail asymptotics from the moment 
asymptotics, but it is natural to make a conjecture. 

Conjecture 3.2. A has a density function /^(x) satisfying 

2 • 3I/6 



fA{x) 



r(2/3) 



00. 



In fact, if A has a density with fA{x) ~ ax^e"'^^'' for some constants 
a, b, c, d, then it is the only possible choice that yields the moment asymp- 
totics (j2.ip . cf. Janson and Louchard [H]. 

Conjecture 13. 2l mav be compared with similar results for several Brownian 
areas in Janson and Louchard 0], see also Janson [3]. Note that in these 
result for Brownian areas, the exponent of x is always an integer (0, 1 or 
2). It is therefore a small surprise that here, the exponent seems to be 1/3, 
corresponding to the power n^^^ in (|2.ip . 
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4. Preliminaries on point processes 

Let © be a measurable space. (In this paper, & is either an interval of 
the real line or the product of two such intervals.) Although a point process 
S will be regarded as a random set {^i} C 6, it is technically convenient 
to formally define it as an integer- valued random measure Yli ■ Hence, 
denotes the number of points that belong to a (measurable) subset 
j4 C (5. Also, X G H is equivalent to H({x}) > 0. See further e.g. Kallenberg 

A Poisson process with intensity d^, where dfi is a measure on S, is a 
point process H such that S(^) has a Poisson distribution with mean fJ,{A) 
for every measurable ACQ, and H(^i), . . . , H(^jfc) are independent for 
every family ^i, . . . , of disjoint measurable sets. Lemma [4. H is a standard 
formula for Laplace functionals, see for instance [6, Lemma 12.2(i)]. 

Lemma 4.1. If E is a Poisson process with intensity d/j, on a set &, and 
/ : S ^ [0, oo) is a measurable function, then 

Eexp{- J2 fiO} = exp{- / (l - e'^^^')) dfi{x)]. 

Lemma [4.21 on the other hand, is more of a digression. The result follows 
from a standard Gamma integral by integration by parts. (The result can 
also be written as 2r(l/2)AV2.) 

Lemma 4.2. If X > 0, then 




5. Proof of Theorem 12.11 

The set {t : B{t) = 0} is a.s. closed and unbounded, so its complement 
{t : B{t) 7^ 0} is an infinite union of finite open intervals, denoted by 
Iv = igu,du), = 1, 2, . . . , in some order. (The intervals cannot be ordered 
by appearance, since there is a.s. an infinite number of them in, say, [0,1]. 
Fortunately, the order does not matter.) The restrictions of B{t) to these 
intervals are called the excursions of B{t). Let be the excursion during 

The local time L{t) is constant during each excursion. Let Ti, be the local 
time during e^, and let ii, := dy — g^ be the length of e^. It is well-known, see 
Revuz and Yor [l3, Chapter XII], that the collection of pairs {(tj^, fj^)}^]^ 
forms a Poisson process in [0, oo) x (0, oo) with intensity 

dA= (27r^3)-i/2 ^^^^^ 

Note also that, a.s., if the excursion e^^ comes before e^^, then Ty^ < t^^. 

Next, consider a Poisson process {Ti}^-^ on [0, oo) with intensity Adt, 
independent of {B{t)}. Assume that the points are ordered with < Ti < 
T2 < ■ ■ ■ . Then Ti, T2 — Ti, ... are i.i.d. Exp(A) random variables with 
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density function Ae . Furthermore, Ti is independent of {B{t)} and thus 
of {A{T)}. It follows from ([O]) that A{Ti) = T^^'^A and consequently 

Ee--^^) = Ee-"^'^'-^ = Ei;{aT^/^) =X T e-^'i^{as^/^) ds. (5.1) 

Jo 

The times T, are called marks, and an excursion is called marked if it con- 
tains at least one of the marks Tj. The marks {Ti} are placed by first con- 
structing {B{t)} and then adding marks according to independent Poisson 
processes with intensities Adt in each excursion. Thus, given the excursions 
{e^}, each excursion e^, is marked with probability 1 — e~^^'' , independently 
of the other excursions. The Poisson process H := {{tu,£u)} defined by the 
excursions can be written as the union H' U H", where 

E' := {{Ty,iy) : Gy IS unmarked}, 
r!' := {{Ty,iy) : is marked}. 

By the general independence properties of Poisson processes, H' and H" are 
independent Poisson processes with intensities 

dA' := e-^' dA = (27r)"i/2r ^/^e"^^ dr d£ (5.2) 

and 

dA" := (1 - dA = {27r)-'/^r'/^ (l - e~^') dr d£, (5.3) 

respectively. In particular, if the lengths are ignored, the local times of the 
marked excursions form a Poisson process H on (0, oo) with intensity 

/ (l-e-^^)dA = Adr, 
Je=o 

where, using Lemma 14.21 

^ PCX) 

A = / (27r)-i/2r (1 - e-^^) di = V2X. (5.4) 
Jo 

Due to the fact that B(Ti) ^ a.s., there exists a unique excursion e,^* 
that contains the first mark Ti, i.e., Ti G 1^*. Let C ■= -^(^i) = '^u* be 
the local time at Ti (and thus during e^*). Since e^* is the first marked 
excursion, its local time C is the first of the points in the Poisson process H 
and hence 

C~Exp(\/2A). (5.5) 

The restriction of B{t) to the interval [0,Ti] consists of all excursions 
with local time Ti, < t^* = Q and the part of e^* on {g^* ,Ti), plus the set 

[0,n]\[jl, = {t<Ti:B{t) = 0} 

V 
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which a.s. has measure and thus may be ignored. Consequently, since 
L{t) = Ty on ly, 

A{Ti) := f ' L{t) dt = [ L{t)dt+ I \{t)dt 

= rJ, + aTi-g,,):=A' + A". 

u:t^<C, 

The sum defined as A! = 'YIiu-tvkC, '''v^u only contains terms for unmarked 
excursions e^. Thus 

A ^ ^ Tu^v 

(r^,^^)6H': r„<C 

Recall that C, is determined by H" (as the smallest r with (r, t) € H" for some 
t) and that r! and S" are independent. Hence, r! and C are independent. It 
follows from Lemma HTl with © = (0, C) x (0,oo) and f{{T,i)) = ari, that 

E(e-"-^' I C) = expj- r (1 - e-""^) dA'(r,£)|. 

Jt=o Je=o ^ 



By (j5?2|) and Lemma g^l 



/ /°^(l-e-°"^)dA'(r,^) 

Jt=0 -'£=0 



r=0 J£=0 

(27r)-V2 /"^ r(e-^^-e-(^+"^)0^"'^'d£dr 
+ ar — Va) dr 



and it follows that 

E(e-"-^' I C) = exp{^/2AC " ^ ((A + "O^^^ - A^/^) }. (5.6) 



Now consider A" = C,{Ti — g^*). Note that Ti — g^* is the location (relative 
to the left endpoint of the excursion) of the first mark in the first marked 
excursion. Since H is a Poisson process with intensity independent of r, 
the location Ti — g^* is independent of the local time C, of the first marked 
excursion. Further, the joint distribution of {iu*,Ti — g^*) has density 

(A)-Ue~^^(27r)-i/2^~3/2 ^^^^^ < y < £ < oo, 
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where the normahzation constant A is given by ()5.4p . Consequently, 



/ e-"^J'(A)-^Ae-^^(27r)-i/2^-3/2 ^^^^ 

/■oo 

= 7r"V2;^l/2 / g-(A+aC)yy-l/2dy 
4=0 

= Ai/2(A + aC)"^/2. (5.7) 

Again, since H' and H" are independent, and A!' are conditionally 
independent given C,. Thus, equation (15. 6p and (15. 7p yield 

E(e-"-4(^i) I C) =E(e-"-^' | C)lE(e~°-^" | C) 

By (j53]) . C has the density \/2Ae"^^, a; > 0, and it follows that 
]Eg-a^(TO ^ A^^"(A + exp{-|^ ((A + axfl'' - A^/^) } dx. 

Finally, the substitution 

provides the slightly simpler formula 



2V2A' 

The result now follows by a comparison with ()5.ip . □ 
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